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Nonstationary Aharonov-Bohm Effect of Current Correlation by Entangled Electrons 

in Quantum Spin Hall Systems 

Wei Chen, 1 ' 2 Z. D. Wang, 2 f] R. Shen, 1 ^ and D. Y. Xing 1 

1 National Laboratory of Solid State Micro structures and Department of Physics, Nanjing University, Nanjing 210093, China 

2 Department of Physics and Center of Theoretical and Computational Physics, 

The University of Hong Kong, Pokfulam Road, Hong Kong, China 

We investigate spin entangled electron pairs transporting separately in the helical edge channels 
of two quantum spin Hall systems. Both dc and ac Aharonov-Bohm (AB) effects of the crossed 
current correlation are studied. In the dc case, it is found that the visibility of the AB oscillation 
of differential noise is equal to entanglement concurrence, and a 7r phase difference indicates singlet 
and triplet entangled states. More remarkably, when an ac flux with frequency fi is introduced, the 
differential noise exhibits nontrivial step structures as the voltage bias V increases, with the step 
edges being located at eV — nfiQ.. The concurrence can be probed directly by measuring the heights 
of these steps. We also show that the relation between differential noise and concurrence is general 
even when spin dephasing exists. 
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As is known, entanglement reflects a kind of nonlocal 
correlation [TJ [SJ. The research on quantum entangle- 
ment has played an important role in quantum informa- 
tion and computation science [3J- Recently, the creation 
and detection of electronic entanglement in solid state 
systems have attracted much interest, for the large-scale 
implementation of quantum information and computa- 
tion schemes [JJ. Moreover, the crossed Andreev reflec- 
tion in mesoscopic s-wave superconductor systems has 
already been confirmed [5] , which is one of the main pro- 
posals for generation of spin entangled electrons in solids 

Though the idea to utilize entanglement in solid sys- 
tems is elegant, a direct experimental evidence is still 
challenging. Several proposals for spin entanglement de- 
tection in solids have been put forward, including the Bell 
inequality tests 018] and the shot noise in a beam splitter 
setup [9]. The former is based on local hidden variable 
theories [5J and the latter utilizes the relation between 
spin entanglement and antisymmetry of electronic wave 
functions [5]. Most recently, we have exploited another 
method [10| , based on a completely different mechanism, 
so called quantum eraser effect [IT] or complementarity 
principle |12j , by noting the inherent relation between en- 
tanglement and quantum interference. It is found that, 
by properly designing the quantum eraser composed by 
the entangled electron pair, the concurrence, i.e. degree 
of entanglement, can be obtained by the amplitude of 
Aharonov-Bohm (AB) oscillation of the current correla- 
tion [TO]. 

In the original proposal of the quantum eraser method, 
stationary AB effect and pure entangled states are con- 
sidered [TU]. In this paper, we generalize our exploration 
to the nonstationary AB effect and mixed states caused 
by spin dephasing. In particular, we focus on a novel 
quantum spin Hall system (QSH) [13] . for the helical 
property of edge states in QSH makes the proposal more 
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FIG. 1. Illustration of the proposed setup. The entangled 
electron pairs are tunneling separately from the entangler (la- 
beled by E) into the helical edge channels of QSHs. The spin 
up (down) channels are sketched as solid (dashed) lines. The 
narrow regions in both sides bring coupling between the upper 
and lower edge channels, and serve as beam splitters. In the 
right side, a magnetic flux $ is enclosed by the loop formed 
by the upper and lower channels. Four terminals (£1,2,-^1,2) 
measure the electronic currents and current correlations. 



realizable. Firstly, the spin and moving direction of heli- 
cal electrons are bounded, which leads to an all-electrical- 
control of electron spin [5] . Secondly, the back scattering 
is forbidden by the time reversal symmetry, so the spin 
relaxation time is extremely long. 

The proposed experimental setup is shown in Fig. [TJ 
where the entangler in the middle region is weakly cou- 
pled to two QSHs in both sides. When a voltage bias is 
applied between the entangler and terminals, entangled 
electron pairs will tunnel into the helical edge channels of 
QSHs. They may be injected into the same QSH, or they 
may be separately injected into different QSHs. However, 
we here only consider the latter process, because the in- 
jection into the same QSH has only a higher order con- 
tribution to the current correlation in the weak coupling 
limit [3j 14]. In each QSH, there is a narrow region with 
the opposite edge channels being coupled, which leads to 
the scattering between edges. It has been demonstrated 
that by tuning the gate voltage properly, the quantum 



point contact structure can serve as an ideal beam split- 
ter without any back scattering [15) . which is beneficial 
to entanglement detection. We here consider that a time- 
dependant flux $(£) is applied in the right QSH of Fig. 
[JJ which leads to the nonstationary AB effect. 

The entangler can be realized by using superconduc- 
tors OE] or quantum dots [H]. Here we adopt the former 
as an example, while the theory itself is actually general, 
which can be extended to other kinds of entanglers. We 
assume that the superconductor is grounded and a bias 
voltage — eV (< the superconducting gap A) is applied to 
all the terminals £1.2,-^1,2, so the Cooper pairs are split 
into edge channels of the QSHs [17] with a tunneling co- 
efficient )(<1. The entangled states may be expressed 
as 0[H] 
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where 0) is the zero temperature ground state of filled 

Fermi sea with a Fermi energy —eV, and a\ a (e) cre- 
ates an electron with energy e and spin a in the 
i side which satisfies the anticommutation relation 
{ai,a{s),aj a i (e')} = Sij5 a ^'5(s — e'). The helicity of the 
edge states indicates the spin-edge correspondence. For 
example, a^ describes an electron moving left along the 
lower edge with its spin polarized up (Fig. [I]). Though 
we here focus on the singlet and triplet entangled states, 
we keep the parameter k tunable to show clearly how the 
physical results changes continuously from the mostly en- 
tangled states (k — 1/2) to product states (k = 0, 1). For 
the pure state in Eq. (flj), the degree of entanglement can 
be described by the concurrence C = 2-\/k(1 — n) [15] , 
Specifically, when C = 1, the "=p" in Eq. (]IJ) correspond 
to the singlet and triplet entangled states, respectively. 

After transmission into QSHs, the electrons get scat- 
tered at the narrow regions as in Fig. [JJ before they 
finally arrive at the terminals. They may change their 
edges during the scattering processes, while the back 
scattering is completely ruled out by tuning the gate volt- 
ages properly [TS]. Thus the electron operators in L\,R\ 
can be written as 



a-Li — t2LO.Lt + tlLO-Ll, 
o.ri — tmaji^ + t 2R a R \^, 



(2) 



where the coefficients t\^.i describe the amplitudes for 
the same-edge and cross-edge transmissions in the i side, 
respectively. The electrons in the other two terminals 
L 2 , R 2 can be described as well, which are omitted here. 
In order to probe entanglement, we calculate the zero 
frequency current correlation between terminals L\,R\, 
which can be expressed by the noise power 
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where the current operator defined as In(t) — (e/h) x 
/ J dEdE'eW-EWal^E^aniE) and the higher order 
term (Ji 1 (i))(jfl 1 (0)) ~ x 4 can be neglected. By utiliz- 
ing Eq. ([2]) and the anticommutation relation of Fermion 
operators, we can obtain the time dependant current cor- 
relation after some algebra 

(I Ll (t)I Rl (0)) = (W f de f de'e^-e'^x 

ll J-eV J-eV 

nT 2L T 2R + (1 - n)T 1L T 1R T C\ZT 2L T 2 rT 1l T 1r cos ip(t) 

(4) 

where the transmission probabilities are Ti 2i = 
\ti,2i\ 2 and the total phase is defined as ip(t) = 
Arg[t2 L t 1L t2 R tin]. From Eq. A), we notice that the 
current correlation can be tuned by the flux in the right 
QSH (Fig. fTl) , and the strongest interference occurs when 
Ti : 2L = Ti,2,r = 1/2, which will be adopted in the fol- 
lowing. 

In order to investigate both dc and ac AB effects, 
a harmonic flux $(£) = $ a sinil£ with amplitude $ a 
is introduced [19 , so the phase has the form ip(t) — 
ipo + 27r($ Q /<f>o) sin fit with a constant term tpo including 
all stationary phase accumulated by transport, scattering 
and also the dc component of the flux, where $ — he/ e 
is the flux unit. 

Given the periodicity of the phase, ip(t + 2iT/Vt) — ip(t), 
the interference term can be expanded as Fourier series 
by utilizing the identity 
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dt [(i Ll (t)i Rl (o)) - (i Ll (*)) (f Rl (o))] , (3) 



where J n is the Bessel function of the first kind. Inserting 
Eq. ([5]) into Eqs. d3l Kb, and performing all the integrals 
over i, e and e', we finally obtain an elegant expression 
for the differential noise (DN) A = dS/deV as 

2 
A = X 2 ^{l TCcosip [j +2J2 J2n0(eV-nm)]}, 

n>0 

(6) 

where we have omitted the variable 2n$f- of the Bessel 
functions for simplicity and 6 is the Heaviside stepfunc- 
tion. Eq. (pi is a main result of the present work, which 
contains important information for probing the concur- 
rence. 

DC case. First we consider the case $ a = when the 
ac component of the flux is absent. In this case, we have 
Jo = 1, J 2n = in Eq. ^, so the expression is just an 
ordinary AB oscillation of the DN as A = x 2 e 2 //z(l =F 
Ccost^o) [10], which is sketched by Fig. [2] for the triplet 
entangled states. The relation between entanglement and 
current correlation is reflected by the equality of the vis- 



ibility defined as v = (A r - 



'^min ) / \^max 1 ^min ) and 



the concurrence, namely v — C. The mechanism of this 
equality is the quantum eraser effect [TT]. Noting that 
the original state by Eq. (IT]) has the typical form ex- 
pressed simply as y/R\ f|) T \/l — ft| |t)i where the right 
electron transports in the interference loop and the left 
electron is the entangled partner. There are only two 
transport channels contributing to the interference with- 
out any spin degeneracy, as a result of the helical prop- 
erty of electrons in QSH. Now if only the right electron 
is measured, there is no AB oscillation of current in Ri. 
The reason is that the left electron serves as the which 
way marker of the right one. For example, if the left 
electron is in the state | f), then we can infer that the 
right electron picks the lower channel. According to the 
complementarity principle |12j , whenever the which way 
information is marked, the interference will not occur. 
However, if we make a correlated measurement between 
the two electrons, or equivalently, the current correlation 
between L\,Ri, then the result will be distinctly differ- 
ent. The fantastic feature of entanglement lies in that as 
long as the spin of the left electron gets measured, the 
right electron will collapse into a corresponding state, 
which may contain no path information at all and the 
interference will recur. To understand this point clearly, 
we can take ti2L — l/v2 as an example, and the spin 
state detected by L 1 terminal is \L\) = 1/a/2(| f) + | 1)), 
with the original spin information | f), | \) totally erased. 
Now if an electron is detected by L\ , we can infer from the 
original entangled state Eq. (fTl) that the amplitude of its 
entangled partner detected by Ri is feiJV^TiiRvl — n, 
where the interference recurs. It is also clear that the 
amplitude of AB oscillation gives out exactly the degree 
of entanglement, because the quantum eraser effect is a 
direct result of nonlocal quantum correlation from entan- 
glement [10] . 

AC case. Then we set (po = to investigate the non- 
stationary AB effect, and the DN has the form A = 
X 2 e 2 /h{lTC[J + 2J2 n>Q J2nO(eV-nhn)]}. The main 
characteristic of the result is the step structures as shown 
in Fig. [3] As long as the the voltage bias V increases 
by hCl/e, the DN jumps onto a new step. The increase 
(decrease) of the nth step is 2C J 2 „ (n > 0) for triplet 
(singlet) state. Remarkably, through the measured ra- 
tios of heights of these steps, one can first obtain the 
flux amplitude $ a . Then with the fixed <& a , the con- 
currence can be measured directly through the heights 
of these steps. In Fig. [3](a), we present the step struc- 
tures of A for triplet entangled states (C = 1), when the 
flux amplitude $ a takes different values. The figures of 
the singlet states are simply the mirror symmetry about 
X 2 e 2 //i, which comes from the "=f" in the original states 
in Eq. M, indicating that they are distinguishable by 
the step shapes. We also note that when the voltage bias 
is big enough, the height of the steps will approach and 
2x 2 e 2 /h for singlet and triplet states, respectively. This 
is not surprising if we notice the identity Xm=-oo ^ n = ^" 




FIG. 2. Plot of the DN as a function of stationary flux tpo 
for triplet entangled states (The singlet states are the mirror 
symmetry about y^e 2 /ft,). The solid lines stand for pure states 
with the dephasing parameter 7=1 and the dashed and 
dotted lines represent the mixed states with 7 = 0.5, 0.3. 
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FIG. 3. Plot of the DN as a function of voltage bias for 
the triplet entangled states (The singlet states are also the 
mirror symmetry about y^e 2 /ft.), (a) The pure triplet en- 
tangled states with the amplitudes of the ac fluxes $ a = 
0.5$o, (fro, 1.5$o corresponding to the lines 1, 2, 3, respec- 
tively, (b) For a fixed amplitude of flux $ a = $0, with differ- 
ent dephasing parameter 7 = 1, 0.5, 0.3. 



Spin dephasing It is worthwhile to discuss the spin de- 
phasing effect in the QSHs. The spin flipped back scat- 
tering is totally suppressed in the helical edge channels 
by the time reversal symmetry, so theoretically, the spin 
relaxation process can be neglected in the present model 
with nano-scale. On the other hand, the spin dephasing 
process always exists, and due to the helical property, the 
spin and edge freedoms are bounded together. Therefor 
any interaction between edge and environment that in- 
troduces a phase uncertainty on the edge freedom will 
also lead to the spin dephasing. Here we would like to 



give a general discussion on the spin dephasing effect and 
its effect on the AB effect of current correlation, without 
limiting in a specific environment, such as a free electro- 
magnetic field and the charge fluctuation in QSH, etc. 

When the electron interacts with the environment, 
they together are generally in an entangled state. The 
interference of the electron can be obtained by per- 
forming integral over all the freedoms of the environ- 
ment. In the quasi-classical approximation, one finds 
that the effect of the interaction is only to multiply the 
interference term in Eq. Q by a scalar product of 
the initial and final states of the environment as |20j 
7 = ^ \fexp[-i/hJ^dtVi(t)]\^ ), during the time r 
of electron traversing through the system, where |£o) is 
the initial state of the environment and Vi(t) is the in- 
teraction term in the interaction picture. Usually the 
phase accumulation in many-body environment satisfies 
Markov approximation, i.e. [Vi (t) , Vi (f )] = 0, for a short 
memory time of environment compared with tq. By us- 
ing the central limit theorem, the dephasing parameter 
is an exponential decay with the phase uncertainty as 
7 (r ) = e" 1 / 2 ^^ 20J. Given that the nonstationary 
case discussed here, it is necessary to compare two time 
scales, i.e. the electron traversing time tq and the period 
of the ac flux 2ir/£l. Only when firo <C 2tt is satisfied, 
the phase <p(t) can be treated quasi-classically as a func- 
tion of time in Eq. Q. Taking the HgTe quantum wells 
as an example [13], the Dirac velocity of edge states is 
given by v ~ 5.5 x 10 5 m/s, so the electron spends about 
t = l/v ~ 1.8ps to transport through the QSH for a 
nano-scale system with a typical length I ~ 1/im. As a re- 
sult, the quasi-classical condition is S7 <^C 3.5THz, which 
is fulfilled by the usual microwave frequency £1 ~ 1GHz 
or by the time-dependant gate voltage on the edge chan- 
nel [H], which results in the ac AB effect. 

Based on the above discussions, the dephasing effect 
can be effectivey described by the dephasing parame- 
ter 7 within the interval < 7 < 1, which changes 
the pure states into mixed states. For the singlet and 
triplet entangled states in Eq. (HI, the original pure 
states l/v2(| f-l) T I it)) 1S now substituted by the den- 
sity matrix p = 1/2 [| U)(U I + I It) (It I T 7(1 tlXIT 
I + I 4-tXtl I)]) witn the concurrence C = 7 [21 . The 
current correlation is now calculated through the for- 

mula (/(t)J(O)) = X 2 i: V eV dej: V eV de'Tv[pi(t)i(0)}, and 
finally we can obtain exactly the same equation Eq. (|6]), 
merely with the concurrence C = 7. Therefor, we con- 
clude that no matter the spin dephasing exists or not, the 
DN formula Eq. ([6| is generally applicable, with C = 1, 7 
corresponding to pure and mixed states. The numerical 
results can be found in Fig. [2] and Fig. J3^b) for dc and ac 
cases, respectively. In the dc case, the visibility v is sup- 
pressed by the spin dephasing effect, while the equality 
v = C between visibility and concurrence is general, for 
both pure and mixed states. In the ac case, the height of 



each step is also suppressed by the dephasing parameter 
7, while all the arguments for pure states are still valid, 
only with C = 1 being replaced by C = 7. Therefor, the 
dephasing effect in QSHs can also be measured through 
the dc and ac AB effect of current correlation. 

In summary, we have studied both dc and ac AB ef- 
fects of the current correlation contributed by entangled 
electron pairs in QSHs. It has been shown that for singlet 
and triplet entangled states, the degree of entanglement 
can be probed by the visibility of the stationary AB os- 
cillation, or the step structures in the nonstationary case, 
no matter the spin dephasing exists or not, providing a 
new method for entanglement detection. 
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